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ABSTRACT 



Context. X-ray surveys of Active Galactic Nuclei (AGN) indicate "cosmic downsizing", with the comoving number density of high- 
luminosity objects peaking at higher redshifts (z ~ 2) than low-luminosity AGN (z < 1). 

Aims. We test whether downsizing is caused by activity shifting towards low-mass black holes accreting at near-Eddington rates, or 
by a change in the average rate of accretion onto supermassive black holes. We estimate the black hole masses and Eddington ratios 
of an X-ray selected sample of AGN in the Chandra Deep Field South at z < 1, probing the epoch where AGN cosmic downsizing 
has been reported. 

Methods. Black hole masses are estimated both from host galaxy stellar masses, which are estimated from fitting to published optical 
and near-infrared photometry, and from near-infrared luminosities, applying established correlations between black hole mass and 
host galaxy properties. Both methods give consistent results. Comparison and calibration of possible redshift-dependent effects is also 
made using published faint host galaxy velocity dispersion measurements. 

Results. The Eddington ratios in our sample span the range ~ 10~ 5 - 1, with median log (Lboi/^Edd) = -2.87, and with typical black 
hole masses M Bii ~ 10 8 M Q . The broad distribution of Eddington ratios is consistent with that expected for AGN samples at low 
and moderate luminosity. We find no evidence that the CDF-S AGN population is dominated by low-mass black holes accreting at 
near-Eddington ratios and the results suggest that diminishing accretion rates onto average-sized black holes are responsible for the 
reported AGN downsizing at redshifts below unity. 

Key words, accretion; galaxies: active 



1. Introduction 

One of the most significant recent discoveries in studying the 
cosmological evolution of AGN has been the discovery that 
lower-luminosity AGN (Z/2-iOkeV ~ - 10 43 ergs~') peak 
in their comoving space density at lower redshift, z $ 1, than 
higher-luminos ity AGN (Zg-iokeV ~ 10 45 - 10 47 erg s~'), which 
peak at z ~ 2 dCowi e et al. 2003; Steffen et al. 2003; Ueda et al 
2003; Bargeretal. 2005; Giacconi et al. 2002; La Franca et al 



2005; Hasinger et al. 2005; Hopkins et al. 2007). 



The implication of the phenomenon for black hole growth 
can be see n from the emissivity of AGN in different luminos- 
ity ranges (Hasinger et al. 2005, Fig. 5b): at redshifts z > 2 the 
major contribution to the total emissivity comes from high lu- 
minosity AGN (L .5-2keV ~ 10 44 - 10 45 ergs" 1 ), while at z ~ 
the major contribution shifts to lower luminosities (Lo.s-2keV ~ 
10 43 - 10 44 erg s~'), with a significant contribution from low lu- 
minosity AGN (Lo.5-2keV ~ 10 42 - 10 43 ergs _1 ). If we simply 
relate the emissivity to the accretion rate onto black holes, it fol- 
lows that a significant contribution to the black hole growth has 
shifted from high luminosity objects at high redshifts to low lu- 
minosity ones at low redshifts. This shift has been named 'AGN 
cosmic downsizing', by analogy with the phenomenon observed 
in the cosmic star forma tion rate density in normal galaxies 
('cosmic downsizi ng', e.g. lCowie et al.lfl996tlBauer ; t al. 2005; 
Pant er et aill2007l) . There, the decline in star formation rate at 
redshifts below z < 2 has been shown to be due to star formation 



preferentially stopping in high mass galaxies towards the present 
day, and the major contribution to the star formation density at 
z ~ arises from low mass systems. The relationship of galaxy 
downsizing to AGN downsizing is unclear: the phenomenon in 
AGN is observed in luminosity, and so an interpretation of the 
phenomenon as being downsizing in the mass of active black 
holes requires either an assumption of the relationship between 
AGN luminosity and black hole mass as a function of redshift 
or measurement of that relationship. Supporting evidence that 
low-luminosity AGN at low redshift are preferentially powered 
by black holes with lower masses than the typical distribution 
of local black holes has been found when using [OIII] emission 
as a measure of activity dHeckman et alj|2~004l) . although when 
radio activity is mea sured instead, hi gher-mass black holes seem 
more active locally dBest et al]|2005l) . 

In this paper we present a determination of the distribu- 
tion of Eddington ratios for an X-ray selected sample of AGN 
in the CDF-S, at z < 1. This redshift range is of particular 
interest, as it is the epoch in which the total AGN emissiv- 
ity is dominated by moderate and low luminosity AGN — the 
AGN downsizing epoch. The Eddington ratios are found by first 
estimating black hole masses from photometrically-estimated 
stellar masses and luminosities. These black holes masses are 
then compared w i th the host galaxy velocity dispersions from 
Ivan der Wei et a fl d2005l) . The Eddington ratio is then calculated 
as A = Lboi/^Edd where Lboi is the AGN bolometric luminos- 
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ity inferred from the observed X-ray luminosity and Lem is 
the Eddington luminosity. We assume Hq = 71 kms~' Mpc" 1 , 
Qm = 0.3 and Qa = 0.7 throughout. 



2. The AGN sample 

The CDF-S has been the subject of a 1 Ms Chandra observation 
reachi ng a flux limit of 4 .5 x 10 16 erg cm~ 2 s _1 in the 2-10 keV 
band (Ro sati et al.l l2002). At redshift z = 1 the flux limit corre- 
sponds to a hard X-ray luminosity of 1.2-iOkeV = 10 424 erg s~\ 
meaning that the AGN population that dominates the downsiz- 
ing, with L 2 -iokeV ~ lO 41 ' 5 - 10 43 - 5 ergs-' at z < 1, is ob- 
served in the surve y. Optical imaging an d identifications have 
been published by Giac coni et alj d2002l) with optical spectro- 
scopic follow-up by ISzokolv et alj d2004l hereafter S04), pro- 
viding spectroscopic redshifts, in addition to X-ray and opti- 
cal classifications. The total number of X-ray sources in the 
catalogue is 347, of which 251 are detected in the hard X-ray 
band. Av ailability of the spectroscopic (S04) an d photometric 
redshifts (IZheng et alJl2004t iMainieri et al.l l2005) together with 
the deep X-ray exposu re has allowed the X-ray spectroscopic 
analysis published by iTozzi et all d2006h . For bright sources, 
these authors fitted the photon spectral index F and the absorbing 
hydrogen column density Nh simultaneously, while for fainter 
sources the photon spectral index was fixed at the mean value 
(r> = 1.8. The Nh distribution was found to have a log-normal 
shape, peaking at log(A r H/cm -2 ) = 23.1. From the fitted spec- 
tra, absorption-corrected X-ray lumin osities were estimated and 
have been made publicly available by lTozzi et all There are 321 
objects classified by their absorption-corrected total X-ray lumi- 
nosity as AGN (Lx > 10 41 erg s -1 ), and the completeness of the 
spectroscopic (X-ray) catalogue with respect to the photometric 
one is ~ 99%. 

Deep BVIz optical photometry in CDF-S is provided by 
the G OODS survey, from ACS/HST imaging dGiavalisco et alj 
120041) with sensitivity limits 27.8, 27.8, 27.1, 26.6 (AB mag- 
nitudes, lOcr limits, in 072 aperture). The deepest R-band 
data avail able (Rab ~ 26.7) are included in the 1 Ms 
catalogue dGiacconi et al.1 120021) . In the near-infrared (NIR), 
lOlsen etalJ d2006l) have re -reduced the ESO Imaging Survey 
data dArnouts et al]|200ll) . with the observations reaching me- 
dian limiting magnitudes Jab ~ 23.1, Kab ~ 22.2 (5<x lim- 
its, in 2" aperture). In the infrared, the Spitzer Wide-area 
InfraRed Extragalactic survey (SWIRE) pr ovides observation s 
in four IRAC bands: 3.6, 4.5, 5.8, 8.0 /mi (tSurace et alj|2005h . 
The 90% completeness limits for SWIRE Data Release 2 are 
m AB = (21.0,20.96,19.84,19.53). CDF-S is also one of the 
fields covered by the COMBO- 17 survey , with photometry in 12 
medium-band and 5 broad-band filters dWolf et alj|2004l) . The 
COMBO-17 magnitude limits are 24.7, 25.0, 24.5, 25.6, 23.0 in 
U, B, V, R, I bands respectively (AB magnitudes, 10cr limits, 
~ 1"5 aperture), while limiting magnitudes in the 12 medium- 
band filters are in the range itiab = 22.0 - 23.9. 

Recently the deepest CDF-S optical to infrared public data 
has been compiled into a uniform photometric catalogue of 
NIR-selected objects by the GOODS- MUSIC (MUltiwave length 
Southern Infrared Catalog) project dGrazian et alj 12006). The 
catalogue is based on imaging in BVIz from ACS/HST, JHK S 
from ISAAC/VLT, U band from the 2.2ESO and VLT-VIMOS, 
and 4 IRAC/Spitzer bands. The authors have developed a 'PSF- 
matching' (point spread function) algorithm, which allows find- 
ing precise object colours between ground-based and space- 
based images. 



In this paper we define and analyse two samples drawn from 
the above data. We limit the redshift to z < 1, as we are inter- 
ested in p robing the epoch f or which AGN downsizing has been 
reported (Barg er et al.l2005l) . Given the X-ray survey flux limits, 
we are guaranteed to probe the AGN population predominantly 
responsible for the inferred downsizing. We include all z < 1 
objects X-ray class ified as AG N (Lx > 10 41 ergs -1 ). There are 
150 objects in the lTozzi et al.1 catalogue that satisfy our redshift 
and X-ray luminosity criteria. 

Objects in our 'spectroscopic' sample are those selected to 
be AGN with secure redshift measurement, secure optical iden- 
tifications and drawn from a catalogue with uniform photom- 
etry. Analysis of the sample uses the BRVIzJHK, and 3.6 /mi 
photometry from the GO ODS-MUSIC catalogue. The catalogue 
has been matched to the ITozzi et ail AGN, with a 4 arcsec ra- 
dius search, and each match has been manually checked. For 
objects with photometric redshifts where identification was not 
secure, the nearest of any multiple matches is kept. The selection 
is restricted to objects with detections in at least 4 bands. The re- 
sulting catalogue comprises 83 objects, from which 55 satisfy 
our 'spectroscopic' criteria, i.e. have secure optical identifica- 
tion and s pectroscopic r eds hift determ ination (redshift quality 
flag 3.0 in Szo kolv et al.l and lTozzi et all catalogues). 

The second sample analysed will contain the remaining ob- 
jects without a secure identification or spectroscopic redshift. 
To obtain coverage of as many objects as possible, it is supple- 
mented with a combined catalogue consistin g of the COMBO- 
17 BVRI b road-band photometry, S ofl JHK, dOlsen et al]|2006l: 
iMov et al.ll2003l) and IRAC 3.6um dSurace et alj2005l) pho tom- 
etry, and the photometric reds hifts fromlZheng et alj (12004). The 
catalogues are matched to the lTozzi et alj AGN. with a 2 arcsec 
search radius, keeping the nearest of any multiple matches. We 
find matches for 41 more objects with detections in at least 4 
bands. The 'photometric' sample now comprises 69 objects. The 
main purpose of creating this 'photometric' sample is to check 
for possible bias that might be caused by restricting the sample 
to objects with secure optical identifications and spectroscopic 
redshifts. This will be described in more detail later. 

In summary, our spectral-template-fitting catalogues have in 
total 124 objects, i.e. 83% of the AGN in CD FS with spec- 
trosco p ic or optical photom etric redshifts z < 1 Szokol v" ;t al.l 
d2004l) : IZheng etalJ d2004l) . Of these, 55 comprise the 'spectro- 
scopic' and 69 the 'photometric' sample. 

Absorption - correc ted X-ray luminosities are adopted from 
the ITozzi et al.l d2006l) catalogue. Our analysis includes four ad- 
ditional AGN from the RDCS 1252-29 field which also have 
published measurement s of host galaxy velocity dispersion 
dvan der Wei et al J 12005) and for which we have obtained new 
X-ray flux measur ements from the Chandra archival images 
dRosati et al.1 120041) . The new fluxes are given in the Table [T| 
and they correspond to absorption-corrected fluxes assuming 
the mean CDF-S intrinsic absorption hydrogen column density 
log(A r n/cm _2 ) = 23.1 and powe r-law spectrum with the photon 
index F = 1.8 dTozzi et al.ll2006l) . Galactic absorption hydrogen 
column density in the direction of the CL 1252 cluster is taken 
tobe JV H = 5.95 x 10 20 cnr 2 . 

The aim is to find Eddington ratios for our AGN sample, 
for which we need the black hole masses and bolometric lumi- 
nosities. Absorption-corrected hard X-ray luminosities together 
with bolometric corrections will give the best available estimate 
of bolometric luminosity. Black hole masses are more difficult 
to estimate: in this paper we infer masses from available host 
galaxy data and use established (at low redshift) correlations of 
supermassive black hole masses and host galaxy properties. The 
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ID 


a 


5 


^2-10keV 




(J2000) 


(J2000) 


(10 -16 erg cm" 2 s" 1 ) 


CL 1252-1 


12 52 45.8899 


-29 29 04.5780 


15.6 ±3.9 


CL 1252-3 


12 52 42.4793 


-29 27 03.5892 


14.2 ±3.6 


CL 1252-5 


12 52 58.5202 


-29 28 39.5256 


5.83 ±2.88 


CL 1252-7 


12 53 03.6396 


-29 27 42.5916 


11.3 ±3.3 



Table 1. Absorption-corrected X-ray fluxes of four CL 1252 
cluster galaxies, measured from Chandra archival images. 
Fluxes were calculated using the galactic absorption hydro- 
gen column density A^h = 5.95 x 10 2() cirT 2 in the direc- 
tion of the CL 1252 cluster, mean CDF-S intrinsic absorption 
logCA^H/cm" 2 ) - 23.1, power-la w spectrum with the photon in- 
dex r = 1.8 dTozzi et alJ l20Q6). and are given in the observed 
2-10 keV band. The optical counterpart coordinates are listed. 
Match of the X-ray sources to the optical is within 2". 



deep multi-wavelength data available for the field allows mea- 
surements of the host galaxy properties that are needed for this 
approach. 

3. Black hole mass estimates 

3. 1 . Host galaxy stellar mass and luminosity 

We have investigated two routes for estimating black hole 
masses for our AGN: (1) finding stellar masses for the galax- 
ies and applyin g a relation between galaxy stellar mass and 
black hole mass (iFerrarese et al.l2006l) : (2) finding the evolution- 
and ^-corrected K-band absolute magnitude and using the well- 
established correlat ion between black hol e mass and K-band 
bulge luminosity dMarconi &Huntl |2003). Both are indirect 
methods, and using two different methods should help demon- 
strate the robustness of the general approach, i.e. finding the su- 
permassive black hole masses from galaxy properties. 

In method (1) galaxy stellar masses are found by fitting 
spectral templates to the measured spectral energy distribu- 
tions (SEDs). SED fitting requires precise relative magnitudes 
(colours), but for inferring galaxy stellar mass the absolute nor- 
malisation of the SED is also required. Furthermore, differ- 
ent catalogues give magnitudes measured in apertures of differ- 
ing diameter, and from images with different seeing. To be as 
close as possible to a un iform catalogue, we use the SExtractor 
dBertin & Arnoutslll996l) automatic aperture (Kron-like) mag- 
nitudes, which provide the most precise total magnitudes for 
galaxies. The flux loss at faint magnitudes is estimated to be 
~ 10% with automatic magnitudes, compared to ~ 70% with 
isophotal or ~ 20% with corrected isophotal magnitudes. We 
use photometry in the BRVIzJHK, and 3.6 pm bands. The longer 
wavelength IRAC bands are discarded because of possible con- 
tamination by dust-reprocessed AGN light. The IRAC 3.6 fim 
band is included, as it will be close to the rest-frame near- 
infrared bands for many of our objects, and NIR light is cru- 
cial for obtaining the correct stellar mass of the galaxy. U pho- 
tometry is not used, again because of possible non-stellar con- 
tinuum contributions. If there is further contamination with the 
rest-frame UV AGN light, we expect to obtain a bad fit to the 
template galaxy spectra: the possibility of a scenario where a 
reasonable fit and non-stellar continuum conspire to indicate a 
younger stellar population is discussed later. 

We use two dif ferent sets of templa t es fr om stellar population 
synthe sis model s : iBruzual & Chariot] d!993l hereafter BC) and 
Mara stonl (120051 hereafter M05), to test the uncertainties com- 
ing from this part of the method. For both sets we allow different 
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types of templates, variation in age of the stellar population and 
reddening. M05 templates include instantaneous bursts (SSP), 
exponentially declining (t = 0.1, 0.3, 1, 2 Gyr), truncated 
(ftrunc =0.1, 0.3, 1, 2 Gyr) and constant star formation histo- 
ries; allowed metallicities are Z Q /5, Z Q /2, Z , 2Z Q ; reddening 
is varied in the range < Ay < 3; initial mass function (IMF) is 
Salpeter (1955). BC templates include instantaneous burst, ex- 
ponent ial (t = 1, 2, 3, 15, 30 Gyr) and constant star formation; 
IMF is M iller & Scald d 19791). metallicitie s are allowed to vary. 
The reddening law is ICalzetti etaf] d2000l) for both sets of tem- 
plates. 

The redshifts are fixed to the spectroscopic values (S04). For 
fitting and stellar mass estimation we us e a modified version of 
the Hyper-Z code (Bolzonella et al. 2000). Both sets of templates 
produce similar stellar mass estimates. 

Black hole mass and the galaxy stellar mass are related us- 
ing the correla t ion i nvolving the dynamical galaxy mass from 
IFerrarese et ai1d2006l) : 

l0g ( W) = (8 - 4? * °-° 8) + ( °- 91 * °- n) l0g (lpfee) ' (1) 

We discuss the validity of using stellar mass as a proxy for the 
dynamical mass in the next section. 

For the second method, we note that finding absolute K- 
band magnitudes at higher redshifts would involve applying k- 
corrections and, more importantly, significant evolutionary cor- 
rections to the measured K-band magnitudes. To make this pro- 
cess less arbitrary, we can utilise the rest-frame absolute K-band 
magnitudes, which are found as a byproduct of the galaxy tem- 
plate SED fitting procedure in the first method, and already con- 
tain the evolutionary and ^-corrections. We estimate the black 
hole masses by ap plying the K-band bu lge luminosity-black 
hole mass relation (Mar coni & Huntir2003l) 

logMBH = (8.21 ±0.07) + (1.13 + 0.12)(logL 0ul - 10.9), (2) 

where the mass and the luminosity are expressed in solar units. 

The galaxy mass estimates from photometry rely on the 
galaxy light not being contaminated by the AGN light. As al- 
ready noted, we do not use the infrared bands beyond 3.6jt/m 
as they are likely to be affected by a detectable contribution of 
dust-reprocessed light. There are further indications that we are 
indeed dealing with the galaxy light only, i.e. that there is no 
significant non-stellar continuum component in the considered 
bands. For 112 (118 for the K-band luminosity method) out of 
124 objects in our input catalogue we are able to obtain a rea- 
sonable fit (in terms of^ 2 ) with a galaxy template, and the shape 
of the probability distribution for stellar mass and K-band mag- 
nitude is such that it allows us to measure their values and errors 
with reasonable confidence. Furthermore, contamination in blue 
bands would make an object appear younger, and have a lower 
stellar mass estimate, compared to a pure galaxy-light object. On 
the other hand, contamination in the near-infrared bands would 
boost the NIR light, making the mass estimate from the K-band 
higher. Thus, we would see a large difference in estimates from 
the stellar mass method and K-magnitude method. 

The fitting procedure does not unambiguously determine the 
morphological class of the object, but it does provide the most 
likely star formation history for the galaxy and the age of the 
stellar population. For more than half of the objects we find that 
the age of the population is larger (often several times) than 
the characteristic timescale for the last star-formation episode. 
About 20% of the objects have the estimated age smaller than 
the characteristic timescale for the galaxy template, and thus ap- 
pear to be star-forming. They do not show a large discrepancy in 
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mass estimate from the two methods, leading us to conclude that 
this is not due to AGN-light contamination. Low or non-existing 
star formation rate in most objects leads us to believe that it is 
reasonable to approximate bulge light with total light, as the in- 
ference is that these are early-type galaxies. When applying the 
relation (O we assume that this approximation holds: the true 
bulge fraction would need to be significantly less to have any 
large effect on our conclusions. In any case, we apply an empiri- 
cal calibration in the next section that should, on average, factor 
out this uncertainty. 



3.2. Comparison with velocity dispersion measurements 

This indirect way of determining black hole masses from corre- 
lations with galaxy stellar properties results in black hole masses 
with uncertainties of the order of the intrinsic scatter in the cor- 
relations. But in addition, systematic offsets could be introduced 
if there is redshift evolution in those relationships, given the rel- 
atively high median redshift for our sample, the indirect way of 
arriving at the final result and the expectation that these relations 
should evolve with cosmic epoch. One of the best-established 
proxies for black hole mass at low redshift is the bulge velocity 
dispersi on. Its possible redshift evolution has received much at- 
tention dWoo et al.ll2006l:lTreu et alj2004 . although not yet with 
definitive results. Here, we take the view that this is the best es- 
tablished relation, with low intrinsic scatter, and unknown red- 
shift evolution, and in this section we calibrate our other two 
black hole mass estimators (Md yn and ^bui) onto the same scale 
of black hole mass at z ~ 1 that is obtained from the velocity 
dispersion. This 're-calibration' is performed using a sample of 
high redshift galaxies with known velocity dispersion measure- 
ments. Then, we can either assume that the Mbh - cr relation is 
assumed not to evolve in redshift range z - — 1 , or we ca n ap- 
ply the evolution in Mbh - cr obtained by I Woo et alj d2006l) and 
iTreu et ail d2004l) . In terms of the results we obtain in section|4] 
it turns out that the former assumption is the more conservative 
(see the later discussion on this point). The shifts in black hole 
mass introduced by this process are relatively small. 

To perform t he re -calibration, we use the sample of 
Ivan der Wei et al.l (|2005), who have made spectroscopic veloc- 
ity dispersion measurements for 25 non-AGN galaxies and 4 
AGN host galaxies in the CDFS field, of which 23 have early- 
type morphology. They also measured 5 non-AGN and 4 AGN 
in the RDCS 1252-29 field. This galaxy /AGN sample covers the 
redshift range 0.62 < z < 1.13, with median redshift 0.97, 
and allows us to verify the applicability of our method and re- 
calibrate the galaxy-black hole relations that we are going to 
use (zero-points only, given the small size of the sample). There 
are at least three reasons for re-calibrating the relations, coming 
from the uncertainties associated with a higher-redshift sample: 
stellar mass estimates ar e less reliable with increasing redshift 
dvan der Wei et al.l 12006): the reports on evolution in relations 
between the galaxy properties and central black hole masses vary 
significantly — from reports of no evolution to redshift evolution 
so strong that it would have to turn off by z ^ 1 to be consistent 
with othe r constraints on AGN and supermassive black hole pop- 
ulations dHopkins et alJl2056t IWoo et alj|2006t Itreu et alj|2004l 
iMcLure et alj|2006t iPeng et al.ll2006l) : finally, we are measuring 
the total galaxy light, but the local relations are given in terms 
of bulge properties. It is important to emphasise that the major- 
ity of the calibration sample consists of galaxies which have no 
known AGN, and the assumption is the properties of the galaxies 
are representative of the AGN hosts. 



We m atch the available p hotometry in the field with the cat- 
alogue of Ivan der Wei et ail and find stellar mass estimates for 
those galaxies, using both previously described methods. 

Galaxies in the sample have effective radii measurements 
and dynamical mass est imates Md yn = aR e cr 2 /G, with a = 5 
(van derWeletalJl2005h . We compare the reported dynamical 
masses and our stellar mass estimates and find excellent agree- 
ment on average (mean offset at the level of a few percent), al- 
though with a large scatter. Thus we proceed with using the stel- 
lar mass estimate as a proxy for dynamical galaxy mass, keeping 
in mind that the results are to be interpreted in an average sense, 
with uncertainties for individual objects being large. 

We can now estimate black hole masses for the 
Ivan der Wei et al.1 galaxies from the measured velocity 
dispersions using the well-established (at zero redshift) re- 
lation between the black hole mass and velocity dispersion 
dFerrarese et al.ll2006l) : 



log 



M, 



BH 



M, 



-■ (8.48 + 0.07) + (4.41 + 0.43) log ( — — - 1 , (3) 

\224kms~ 1 / 



with a rms of ~0.3 dex in log Mbh- These black hole mass 
estimates may then be compared with our mass estimates us- 
ing the same method as for our AGN hosts, namely applying 
relations ([TJ and to stellar mass and rest-frame, evolution- 
corrected K-band luminosities. This comparison tests for any 
differential evolution between the black hole mass estimators 
and also allows us to place our black hole mass estimates onto 
a common scale defined by the zero-redshift Mbh - cr relation, 
but does not of course take into account any evolution in that re- 
lation, which will be discussed later. Fig. \T\ shows the results of 
this comparison. The differences between the two stellar popu- 
lation synthesis models are smaller than the uncertainties within 
one model, and we plot results from M05 models only. Top 
and bottom panels show the stellar mass method and K-band 
luminosity method results, respectively. K-band luminosity er- 
rors are smaller than stellar mass errors because the rest-frame 
K-band luminosity is well constrained by the 3.6 /mi measure- 
ments, whereas the variation in the total mass yielded by dif- 
ferent templates can be significant. Errors on plotted quantities 
are estimated by identifying a change Ax 2 = 1 in the fits. Dotted 
lines are the nominal relations (Q~|[2|i, and the solid lines are mod- 
ified relations needed to obtain agreement between our model 
predictions and measured velocity dispersions on average. The 
zero-point shifts in the Mbh - M» and Mbh - Lk relations cor- 
respond to black hole mass shifts A log Mbh = 0.09 and 0.31 for 
the stellar mass and K-band luminosity methods, respectively. 
We again assume that the host galaxy velocity dispersion is 
close to that of th e bul ge velocity disper sion, as also assumed by 
IWoo et al.l d2006l) and iTreu et all d2004l) . justified by the obser- 
vation that velocity dispersion depend s only weakly on mea sure- 
ment aperture for early-type galaxie s ( Jorgensen et al.|[l~995l) . As 
already noted by Ivan der Wei et al .1 (2006), we find that stellar 
population synthesis models typically underestimate the velocity 
dispersi ons somewhat. This topic is beyond the scope of this pa- 
per (s ee lvan der Wei et a l. 2006; Drorv et al.ll2004l:lRettura etail 
2006), and for our purposes we simply rescale our mass esti- 
mates to obtain, on average, a match to the measured velocity 
dispersions. 

As noted above, we have thus far assumed that locally mea- 
sured relations (Q~|[2][3|l hold at higher redshifts, and have simply 
shifted the black hole masses by A log Mbh = 0.09 and 0.31 
for the stellar mass an d K-magnitude methods , respectively to 
obtain agreement with Ivan der Wei et al.l (|2005) measurements. 
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Fig. 1. Calibration of our two black hole mass estimate meth- 
ods using the measured velocity dispersions in CDF-S field, for 
a sample of mainly normal galaxies (with 4 AGN in the sample), 
a) Stellar mass method: black hole masses are calculated from 
the measured velocity dispersions using relation @ and com- 
pared to stellar mass estimates, b) K-band magnitude method: 
black hole masses are again calculated from the measured ve- 
locity dispersions using relation (0 and compared to absolute 
K-band magnitudes. Dotted lines are the nominal zero-redshift 
relations (Q~|[2]i. Solid lines are shifted by the amount required to 
make the black hole masses from measured velocity dispersions 
and model quantities agree on average. 



This way, we have factorised out any redshift-related uncertain- 
ties in the stellar light from the host galaxies. Given that the shift 
that has been applied is relatively small, the difference in median 
red shift between the A GN 'spectroscopic' (median redshift 0.7) 
and Ivan der Wei et al.l (median redshift 0.97) samples should be 
of little significance. 

There is however experimental evidence that the black hole 
mass-velocity dispersion relation (0 itself might have redshift 
evolution, in the sense that at higher redshift, galaxies with the 
same velocity dispersion host more mass i ve black ho l es than lo- 
cally (IWoo et al.H2006l iTreu etalJl2004T) . IWoo et al.l (l2006h re- 
port a shift in the black hole mass by A lo g Mbh - .57 ± 0.11 
at z = 0.36, similar to that repo rted by iTreu et all d2004l) . At 
high redshift. IShie Ids et alj (120031) have repo r ted no evolution in 
the Mbh - 0" relation whereas S hields et alj d2006l) do find evi- 
dence for evolution. If evolution does occur, the bootstrapping to 
the zero-redshift Mbh - 0" relationship that we have carried out 
would result in black hole masses being underestimated, and the 
Eddington ratios that we estimate in the following section to be 
overestimated. 




Log [L 2 _ 10 keV /erg s ] 



Fig. 2. Black hole masses and hard X-ray luminosities of the 
'spectroscopic' sample of AGN in the CDF-S. Blue open circles: 
black hole masses estimated using the stellar mass method; red 
stars: black hole masses calcul ated from measured host veloc- 
ity dispersions for the AGN in Ivan der Wei et al.l (120051) . Solid 
squares mark objects which show presence of broad lines in op- 
tical spectra. Lines (dotted and dashed) denote lines of constant 
Eddington ratio of 0.01 and 1 . 0, respectively. R edshift evolution 
in the M BH - o~ relation (e.g. IWoo et ail 1200 6). would cause a 
shift in all points to higher mass by A log Mbh ~ 0.57. 



4. Eddington ratios in CDF-S 

4.1. Analysis of the 'spectroscopic' sample 

Having estimated the black hole masses for the sample, and 
given the measured hard X-ray luminosities, we can now es- 
timate Eddington ratios (the ratio of an object's total luminos- 
ity and its Eddington luminosity corresponding to its black hole 
mass, A = L boi/^Edd)- Bolometric luminosities are obtained by 
applying the iMarconi et al.l (|2004) bolometric correction to the 
absorption-corrected hard X-ray luminosities. Black hole mass 
estimates are described in the previous section. In summary, they 
are calculated using the black hole-galaxy correlations, but with 
zero-points shifted to obtain agreement with the available high- 
redshift velocity dispersion measurements or dynamical mass 
estimates. The absorption-corrected hard X-ray luminosities and 
black hole mass estimates for our sample are shown in Fig. [2] 
together with lines of constant Eddington ratio of 0.01 (dotted 
line) and 1.0 (dashed line). 

Also shown on Fig. [2] are black ho le masses estimated 
directly from the Ivan der Wei et alj d2005) velocity dispersion 
measurements, assuming the zero-redshift relation (fj), for the 8 
AGN with such measurements. It can be seen that these yield 
a similar distribution to the stellar mass estimates, as expected 
from the previous section. This provides independent evidence 
that systematic contamination from AGN light in the photomet- 
ric mass estimation is not a significant cause of bias. 

Finally, we combine our bolometric luminosities and black 
hole masses into Eddington ratio estimates. Fig. [3] (shaded his- 
tograms) shows the distribution of Eddington ratios for our sam- 
ple, using both stellar mass and K-magnitude methods. Median 
(mean) values for the two methods are log A = -2.76 (-2.87) and 
-2.64 (-2.58), respectively. 
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Fig. 3. Distribution of Eddington ratios for the 'spectroscopic' 
sample (shaded histogram) using (a) the stellar mass method, (b) 
the K-magnitude method. Median Eddington ratios are log A = 
-2.87 and -2.58, respectively and the corresponding values are 
marked with arrows. The solid line histogram shows the dis- 
tribution of Eddington ratios for the combined 'spectroscopic' 
and 'photometric' sample. Median Eddington ratios for the com- 
bined distributions are log A = -2.07 and -2.16 (marked with 
arrows). Given the additional uncertainties in the 'photometric' 
sample, the joint sample is expected to be broader than the distri- 
bution for the 'spectroscopic' sample only. The conclusions are 
largely unaffected by the choice of sample. 



4.2. Analysis of the 'photometric' sample 

Initially we have a hard X-ray selected sample, but the reduction 
to the secure redshift and optical-X-ray identification subsample 
could potentially be biasing us to brighter, and thus likely more 
massive objects (as they are more likely to have spectroscopic 
redshift measurement). To test for possible biases of our 'spec- 
troscopic' subsample, we use the photometric redshift estimates 
in CDF-S (IZheng et al.ll2004l) . which are available for ~ 99% of 
the X-ray detected sources, with the average redshift accuracy 
~ 8%. We supplement the non-secure sample with a combina- 
tion of photometric catalogues in the field, resulting in photo- 
metric data for ~ 83% of the whole z < 1 X-ray selected sample. 
The combined sample colours are less precise than those from 
the uniform GOODS-MUSIC catalogue. In addition, neither op- 
tical identifications nor the photometric redshifts are certain for 
all objects, so there are likely to be several objects with erro- 
neous Eddington ratio estimates. This makes the catalogue less 
reliable compared to the 'spectroscopic' one, but still the best 
estimate we can make for the objects missing in our 'spectro- 
scopic' sample. We proceed to find stellar mass estimates for the 
objects. Fits with a well-defined minimum in x 1 are found for 
61 out of 66 objects. For broad-line AGN (18 and 11 objects in 
'photometric' and 'spectroscopic' samples, respectively) stellar 
masses are more likely to be incorrect as the photometry could 
have a significant contribution from the central object, but we do 
not remove them from the sample. 

We find that the distribution of Eddington ratios for all ob- 
jects (solid line histogram in Fig. [3), which includes the 'spec- 
troscopic' sample (shaded histogram in Fig. [3j and the 'photo- 
metric' sample, follows closely the distribution of the 'spectro- 
scopic' sample, convincing us that our subsample is not biased 
towards optically bright objects with high black hole masses and 
low Eddington ratios. Fig. |4] compares the luminosity and red- 
shift distributions of the 'spectroscopic' and 'photometric' sam- 
ples. Slightly higher redshifts (mean Az ~ 0.1) make the 'pho- 
tometric' sample fainter than the optical spectroscopy limits but 
those galaxies are not systematically less massive. 



39 40 41 42 43 44 45 46 0.2 0.4 0.6 0.8 

Log(L 2 _ 10keV /erg s~') z 

Fig. 4. Distributions of absorption-corrected hard X-ray lumi- 
nosities (left panel) and redshifts (right panel) for the 'spectro- 
scopic' sample (shaded histogram) and for the combined 'spec- 
troscopic' and 'photometric' sample (open histogram). There is 
a shift to lower X-ray luminosities and to higher redshifts for the 
'photometric' sample. 



5. Discussion 

5.1. Distributions of Eddington ratio 

The analysis presented here indicates that X-ray selected AGN 
at the 'cosmic downsizing epoch' in fact show a wide range in 
Eddington ratio, with a median value log A ~ -2.8. We now 
discuss how this result compares with other determinations of 
the distribution of AGN Eddington ratios in the literature. 

At low redshift, iPanessa et alj (I2006I) find a similarly broad 
range of Eddington ratio for local Seyfert galaxies selected on 
the basis of optical spectra. They measured nuclear X-ray fluxes 
for the galaxies, and find an active nucleus in all but 4 out of 47 
Seyfert galaxies. The sample covers the bolometric luminosity 
range Lboi ~ 10 41 - 10 44 erg s _1 , very similar to the one for our 
sample, but extending to lower luminosities because of the abil- 
ity to resolve the nucleus in X-rays. The black hole mass range 
(estimated mostly from the Mbh - cr relation) is almost identical 
to our sample at z ~ 0.7. The estimated Eddington ratios for this 
sample again have a broad distribution, in the range log A — -1 
to -7, and median at A ~ 0.01 (type 1 AGN) or ~ 10 3 (all 
AGN). 

However , for br ight (L bol ~ 10 44 - 10 47 erg s 1 ) local AGN, 
I Woo & Urrvl (12002) find Eddington ratios typically in a narrower 
range 0.001-0.1, while at z ~ 0.7 the ratios are somewhat higher, 
0.01-1. For optically selected AGN at 0.3 < z < 4, with bolo- 
metric luminosities in the range Lboi ~ 10 45 - 10 47 erg s _1 , the 
Eddington ratios have a narr ower distribution (0.1- 1), and no 
apparent redshift dependence dKollmeier et al.l |2006). The most 
obvious difference with these samples is the different bolometric 
luminosity range, even at z < 1 . 

The variety of reports on Eddington ratio mean value and 
width of distribution can be understood from the following sim- 
ple illustration of selection effects. Fig. [5] shows the distribu- 
tion of Eddington ratio we would expect to observe for samples 
selected at three differing luminosities. In this illustration, the 
black hole mass function has been tak en to be the double po wer- 
law form for 'active' black holes of iGreene & Ho (2007), the 
intrinsic distribution of Eddington ratios is taken to be uniform 
in log A over the range -5 < log A < 0, and expected distri- 
butions for three choices of bolometric luminosity are shown: 
L = 0.001L*, 0.01L*, 0.1L*, where L* is the Eddington lumi- 
nosity for blac k holes at the 'brea k' in the double power-law 
mass function (IGreene & Holl2007l give the black hole mass at 
the 'break' as M* H = 10 7 32 M o at z = 0). 
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Fig. 5. Expected distributions of Eddington ratio A in the sim- 
ple illustrative model described in the text. Three 'samples' are 
assumed drawn from the active black hole mass function of 
iGreene & Hoi (2007), with a uniform intrinsic distribution of 
log(/l), at luminosities L = 0.00 1L* (red curve), 0.0 1L* (green) 
and O IL* (blue). 



The illustration demonstrates that, while drawn from 
the same (unremarkable) intrinsic distribution, the observed 
Eddington ratio distribution is expected to be strongly dependent 
on the luminosities probed by an AGN sample: selecting low lu- 
minosity objects results in a wide Eddington ratio distribution 
with a low mean value, whereas selecting high luminosity ob- 
jects results in a narrow distribution with a high mean Eddington 
ratio. 

Additionally, the observed distribution of log A is a function 
not only of the form of the intrinsic distribution but also of the 
shape of the black hole mass function and the luminosity selec- 
tion and range covered by the sample. Further selection effects 
may also arise, if there are system atic varia t ions in AGN SEDs 
as a function of A, as indicated by iBorosonl d2002l) and implied 
by the difference in Ed dington ratios typic ally seen for type 1 
versus type 2 AGN (e.g. lPanessa et a l. 2006). 

5.2. Cosmic downsizing 

Similar considerations to those above lead us to conclude that it 
is necessary to understand the distribution of black hole masses 
in AGN samples before physical interpretations may be made 
of the phenomenon of 'cosmic downsizing'. The most striking 
result from the analysis presented here is that the AGN respon- 
sible for the peak in space density at z ^ 1 at moderate AGN 
luminosities cannot be described as being low-mass black holes 
accreting at high rates. This argues against the simplest interpre- 
tation of 'cosmic downsizing': the phenomenon is not due to an 
increasing dominance of low-mass black holes, the typical black 
hole mass in the CDF-S AGN sample is ~ 10 8 M . 

Downsizing thus appears to have a rather more complex ori- 
gin. The typically low Eddin gton ratios we find cou ld be consis- 
tent with the model of lHopkins & Hernquistl d2006l) . where low- 
level AGN activity is fuelled by stochastic accretion of cold gas, 
and dominates the AGN population from z — up to redshifts 
Z ~ 1, and at higher redshifts merger-driven AGN fuelling might 
be dominant. This could be a generic expectation of hierarchical 
structure formation: merger-driven bright phases of AGN activ- 
ity exist in high redshift universe, possibly alongside the less 



bright activity powered by stochastic accretion. At low redshifts, 
once the mergers become less frequent, stochastic accretion be- 
comes the dominant mode. 

More generally, however, if black hole growth is coeval with 
galaxy growth, then we expect the mean accretion rate onto 
galaxies, and hence their assoc iated black holes, to decrease with 
cosmic time dMiller et alj2006l) . Indeed, decrease in mean accre- 
tion rate must occur, since the integrated AGN luminosity den- 
sity decreases with cosmic time at z ^ 2, yet we do not expect 
the integrated mass in black holes to decrease substantially with 
time (rather, the latter should increase as black holes continue 
to grow). Our results suggest that the AGN cosmic downsizing 
at z ^ 1 is not a symptom of 'anti-hierarchical' behaviour, but 
in fact may be a reflection of the process of the dying-down of 
cosmic accretion and a shift in the typical luminosity of mas- 
sive black holes to lower values. Whether this tren d continues to 
still lo wer redshift is still an open question: the iPanessa et alj 
(2006) black hole mass and Eddington ratio distribution ap- 
pear similar to the CDF-S AGN of comparable luminosity, yet 
it seems that lower-luminosity type 2 AGN at low redshift are 
domin ated by lower mass bl ack holes with moderate Eddington 
ratios dHeckman et al.ll2004l) . 

Finally, we note that for a sample of high-redshift (median 
z = 2. 2) sub-mm-selecte d galaxies in the Chandra Deep Field 
North, iBorvs et all d2005l) have found that, if their black holes 
are assumed to be radiating at the Eddington limit, the black 
hole masses in those high-redshift star-forming galaxies are 1- 
2 orders of magnitude smaller than in galaxies of comparable 
mass in the local universe. Given the high redshift of the sample, 
and the discussion of evolution in mean accretion rate above, 
it may well be that typical Eddington ratios are higher than for 
our sample. But interestingly, if we assume that the effect they 
see is actually from the sub-Eddington accretion rates similar to 
the ones we report here (~ O.OlLEddX this would bring the stel- 
lar mass-black hole mass relation for high-redshift star-forming 
galaxies into agreement with the local one. Reliable black hole 
mass estimates at higher redshifts would be a key advance in 
understanding this area. 

6. Conclusions 

We have estimated Eddington ratios for two samples of hard 
X-ray selected AGN in CDF-S with median redshift z = 0.7. 
The primary 'spectroscopic' sample has secure redshifts and op- 
tical identifications and spans the bolometric luminosity range 
Lbo\ ~ 10 41 - 10 45 erg s . The majority of the sources are radi- 
ating at low Eddington ratios in the range A ~ 1 - 10~ 5 , with the 
median log A = -2.87. A larger sample, based on photometric 
redshifts, has fewer selection effects, but larger uncertainties re- 
lated to optical identification and Eddington ratio estimates: this 
fainter sample has Eddington ratios that span the same range, 
and the median for the whole sample is log A = -2.07. Black 
hole masses are in the range 10 5 - 10 9 M , with the distribution 
peaking at ~ 10 s M . It is likely that fainter X-ray flux limits 
would reveal even more sources radiating at low Eddington ra- 
tios. 

We have discussed how the broad distribution of Eddington 
ratios arises because of the relatively low luminosities probed 
by the sample, and that in general observed distributions are 
strongly dependent on the selected luminosity range of AGN 
samples. Based on the estimated Eddington ratios and black hole 
masses for the CDF-S AGN, we argue that diminishin g accretion 
rates onto average-mass supermassive black holes (Mille r~e"t alj 
2006) are the underlying cause of the observed cosmic AGN 
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downsizing at z ~ 0.7, contrary to an interpretation in which 
most of the activity occurs in rapidly-growing low-mass black 
holes. 
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